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Abstract 

We show how to calculate the quantum mass correction to (1+1)D solitonic field theories using 
numerical methods. This is essential if we want to find the corrections to non-integrable models. We 
start with a review of the standard derivation of the first order quantum correction. Then, we re- 
derive a trace formula which allows us to compute the mass correction mode by mode. Specifically, 
we are interested in the extent to which the lowest modes from both, the soliton and the vacuum, 
sectors give the leading contribution. We apply the technique to both the Sine-Gordon and the 
4 -kink model. Then, we compute all the modes numerically and hence the first order quantum 
contribution to the mass of the Sine-Gordon and tfi 4 soliton. 

1 Introduction 

We focus on the computation of the mass correction of (1+1)D solitonic theories theoretically and 
numerically. We follow Rajaraman's semi-classical quantisation procedure; another option is via 
path integrals - see Dashen et al. Q. Specifically, we are interested in the use of the normal modes 
to compute the mass correction and the extent to which the lowest modes from both, the soliton 
and the vacuum, sectors give the leading contribution. We start with a review of the derivation 
of the first order quantum correction. Then, we derive the trace formula from first principles]^] 
We use this formula to re-compute the mass correction and show that the lowest modes are the 
most important ones. Then, we calculate the lowest modes numerically and hence the leading 
contribution to the mass correction. 



2 Mass Quantum Correction: General Idea 

We start out with a general lagrangian for a (l+l)-dimensional field theory 

L=\tf-\cp l2 -V{cp) (1) 

for the scalar field <p{t, x) and with the potential V being positive. The time-independent Euler- 
Lagrange equation leads to 

-0" + ^ = O. (2) 



1 this is not done in 
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We quantise around the minimal-energy static solution (f> s t{x) satisfying (j^), which could be the 
vacuum or the minimal-energy solution in a non-zero topological sector. The semi-classical ex- 
pansion states that the quantum field <p(t,x) is the classical static field <fi s t{x) plus a quantum 
correction field e(t, x), 

4>(t,x) = (f> st (x) + e(t,x) (3) 

where " reminds us to treat the function as a quantum object, an operator satisfying, possibly 
non-commuting, commutation relations with other operators. We have to substitute ([}]) into the 
hamiltonian 

H(<f,)= r dx(U 2 + yi+v^)) (4) 



and obtain, in orders of e, 




= (2) Quantum 

where tt is the conjugate momentum of e. We have used integration by parts, set the boundary 
terms to zero and Taylor expanded the potential term in powers of e. We split our hamiltonian 
into three parts, 

H -^Classical ^Quantum ^-HigherOrder ) (6) 

where the classical mass/energy can be found by substituting <p st into (^). We concentrate on the 
lowest order quantum energy HQ uantU m- This is justified in the framework of a perturbation theory 
where higher order terms are neglected due to their smallness in comparison to the lowest term. 
This approximation is justified if the potential is roughly harmonic around the static solution and 
e™ terms depend on the coupling constant A in the form A™, for example. If these conditions are 
not fulfilled or a better accuracy is wanted, one can resort to well-known perturbation techniques of 
standard quantum mechanics. However, let us emphasise that this would considerably complicate 
our task. Therefore, Barnes et al.'s attempts are already significantly hampered for those multi- 
skyrmions whose potential does not have a steep valley i.e. arc not suitable for a harmonic oscillator 
approximation. This is the case for B = 2, for example.^] 

Our quantum hamiltonian, in the harmonic approximation, has the form 

/ £. d 2y 

2H{t,x) = it 2 {t,x) + e{t 1 x)\-^+- nj )e(t,x) (7) 




where A 2 is an operator. If A 2 acts as a number, the hamiltonian has the form of an harmonic 
oscillator and we know the quantisation procedure. In effect, we have to solve the eigenvalue 
equation 

A 2 e(t,x) = LU 2 e(t,x) (8) 

and this is equivalent to the time-independent Schrodinger equation. We have to decompose the 
quantum field e in terms of a complete set of real and orthonormal cigcnfunctions. Therefore, 

oo 

= ^2 Qn(t)fjn(x) (9) 
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with 



^2vn(x)fj n (y) = S(x-y) 

n 

dx[f) n (x)r) m (x)\ = <5„ im (10) 



and 



A 2 ^) =Uifk(x)- (11) 
We substitute the decomposition (^) into (Q), integrate over x and use the constraints ( |l0| ) on the 
eigenf unctions. We are left with 

2H(t)=J2[pl(t)+^€(t)] (12) 

n 

where p n is the conjugate momentum of q n . The hamiltonian is an infinite sum of harmonic 
oscillators of frequency oj n . We have reduced our quantum field theory with an infinite number of 
degrees of freedom to a pseudo particle quantum mechanics with an infinite number of harmonic 
oscillators. We are now able to use standard particle quantum mechanics. Using the operator 
method with Heisenberg's commutation relation 

[PnAn] = ih, (13) 

we get 



H = hy(a n + l)Lu n (14) 



where a n is the ath energy level of the nth oscillator. Naively speaking, we have all the information 
we need to calculate the mass correction: the classical mass and the quantum correction to first 
order i.e. the quantum hamiltonian. Unfortunately, if we were to calculate the quantum mass in a 
specific model, we would quickly realise that the mass is divergent; the infinite number of oscillators, 
an inherent feature of any quantum field theory, being the cause of this divergent result. In the 
next section, we describe, using the tfi 4 model as our example, how to extract a meaningful quantum 
properties from our naive expression. 



3 Mass Quantum Correction: Derivation 

We follow Rajaraman's procedure. For a very detailed and lucid description, we refer to Rajara- 
man's book j|, section 5.3]. We have re-done all the calculations by hand and with the use of Maple. 
Further we have expanded on the discussion of some parts, given more details and corrected some 
typographical errors.^ 

The hamiltonian of the 4 kink model is 

2H = Jdx (> + r - m V + ^ + |0 (15) 

where m is the mass of the field <j> and A the self-coupling constant. In topological charge sector 
zero, the minimal-energy solution i.e. the vacuum is 

ct> st {x) = ±^= (16) 
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and probably introduced others. 
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and, in topological charge sector one, the minimal-energy solution i.e. the 1-kink is 

mix — a) 



4>st{x) — i^^tanh 
V A 



V2 



(17) 



where a indicates a translational invariance (this will lead to a zero mode). We quantise around 
the solutions with the positive sign in front which is the soliton - compared with the negative sign 
in front which is the anti-soliton. 

The corresponding eigenvalue equation for the vacuum is 



dx 2 



The eigenfunctions and eigenvalues are 



2m 2 



r)i(x) = w^j i)i{x) 



f)k(x) — exp(ikx) 



and 



J V,k 



k 2 + 2m 2 . 



We use periodic boundary conditions in a box of length L and 

k n L = 2im 



(18) 

(19) 
(20) 

(21) 



where n is an integer. The continuum limit is reached by taking L to infinity and any discrete sum 
over k n , or simply n, turns into an integral over k of the form 



E 



dn = 



L 
2^ 



(22) 



using the constraint (|2|) on k n . 

The corresponding eigenvalue equation for the 1-kink is less trivial 



-\& 2+{M2z - 1) 



m~ 



(23) 



where z = mx /y/2 for convenience. In fact, it belongs to a class of special Schrodinger equations; 
the Sine-Gordon model being another example.^] There are two discrete modes; one zero mode 



fjo(z) = cosh 2 z with w^-p = 

and a second discrete mode 

fji (z) = sinh z cosh -2 z with ld 2 k = - 
The continuous eigenmodes, which we label with q, are ^ 



(3 tanh 2 



1 — q — 3iq tanh z) 



with 



K,q 



m 2 ^ + 2 



(24) 
(25) 

(26) 
(27) 



4 private communication from Jackiw and see pages 683-684] 
5 the eigenmodes are Jacobi polynomials in tanh z. M 
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Imposing periodic boundary conditions becomes more tricky. For z — > ±00, we can use the asymp- 
totic form of f] q {z) 



%{z) -> e iqz (2 - g 2 =F 3iq) 

— ► vV+4) 2 (g 2 + l)exp 



(28) 



where S(q) is the phase shift of the scattering states from the viewpoint of the Schrodinger equation 

3q 



S(q) = —2 arctan 



The condition imposed by the periodic boundary is 



2-q 2 



(29) 



q n (^j + 5(q n ) = 2nn (30) 
and the sum over n becomes an integral over q in the limit where L goes to infinity: 

I* + (3!) 

using the constraint ( |30| ) on n. Let us briefly note that the zero mode, which we have found, is noth- 
ing else but the manifestation that the 1-kink solution is translationally invariant. We can imagine 
the infinite-dimensional potential space at its minimal-energy location having a 1-dimensional val- 
ley along which we can move our solution by varying a without changing the energy of the system. 
This zero mode certainly needs to worry us, for our harmonic oscillator approximation assumes 
steep valleys in all dimensions. However, this is only problematic if the zero mode is coupled with 
another mode and this does not happen in the computation of the mass correction up to first order. 
A proper treatment of zero modes is done with collective coordinates: see Rajaraman chapter 
8]- 

We have all the necessary information on the eigenvalues and should naively be able to compute 
the mass of the 1-kink up to first order quantum corrections. Using (0) and the eigenvalues of the 
1-kink solution, we get an expression for the energy 



272m 3 _ V3 U f^~o 




Discrete 



Continuous 



Finite 



Divergent 



which includes the finite classical energy, no contribution from the zero mode due to its zero 
frequency, a finite contribution from the second discrete mode and a sum over the continuous 



modes. Unfortunately, if we were to perform the integral over q using (30), we would find it to be 
divergent. This clearly shows that our naive treatment of quantum field theory is inadequate. To 
have a finite answer, there are two modifications we have to make. 

3.1 Energy level difference 

Let us write out the expression for the vacuum energy up to first order quantum corrections 

Ev = lhJ2Vk n + 2m 2 . (33) 
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Using (|22|), we get 

hL 



E v = — dkyjk 2 +2m 2 (34) 

47r J-oo 

which is a quadratically divergent integral. Thus, even the leading quantum contribution to the 
classical vacuum is not finite. However, we can follow the example of newtonian gravity which 
defines potential energy as the difference between two states. It makes physical sense to define our 
naive vacuum energy, even though it is infinite, as the lowest of possible energy states in the theory 
i.e. to put it equal to zero and hence to subtract it from our naive calculation of the 1-kink energy. 
We get 



E K - E V = E K = + hm^= > J | m\j £g» + 2 - y/k* + 2m 2 ] (35) 

E Sini te 



and we label all finite terms collectively E finite- We go to the continuum limit and perform the 
integral over k. Therefore, we re-express q n in terms of k n 

q n = — k„ — (36) 

m \ LI 



using the boundary conditions ( |30| ) and (21). The expression in the sum takes the form 



rn 



\ql + 1- yjkl + 2m 2 = ^J(k n -^^j +2m 2 - V^ + 2m 2 

0(1/L 2 ) 



k n 5(kn) v 



L^/kl + 2m 2 

~ arctan ( ^_^_") - K + 0(1/ L 2 ) (37) 
L \V2m 2 -kl) jk* +2m 2 



where we have Taylor expanded the first line, used expression (|29|) and Taylor expanded it. Both 
Taylor expansions are in and make sense, for we take the box size L to infinity later. Using (|22]), 
the expression for the energy becomes 

~ 7& f ( 3ttt< k \ k 

E K = E finite + — dk arctan + 0(1/L) . (38) 

2-k J \ V2 m 2 - fc 2 / Vfc 2 + 2m 2 

The dependence on the box size goes away for L — > oo and we are allowed to neglect the 0(1/ L) 
and higher terms. As our notation indicates, we perform integration by parts. The boundary term 
has the form 



— lim 

2"7T a^oo 



arctan | ^ — „ - - „ | \J k 2 + 2m 2 
\^2 m 2 - k 2 J 



(39) 



This limit is ill-defined for trivial substitution of a — oo: the arctan function gives us and the 
polynomial function oo. Therefore, we use l'Hoptial's rule and obtain a finite answer 

Hm -i= (40) 
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which we include in our E 'finite- The integral obtained by integration by parts has the form 

3V ~ 2m fdk - P + m2 . (41) 

2tt J VA: 2 + 2m 2 (2k 2 + m 2 ) V ' 

We put a cut-off A on the k limits and change to the variable p = k/m. We get 

lim p - p2 + 1 (42) 

^°°J± v /^+2 (2p 2 + 1) 

If we perform this integral, we still find a logarithmic divergence plus a finite contribution. We 
need to cancel the divergence with another term. We need to look closer at the infinities produced 
by the infinite degrees of freedom of a field theory. 

3.2 Normal-Ordering and Counter-terms 

We have to normal-order the hamiltonian and introduce counter-terms. We do not give a full 
introduction to all these more complicated ideas and refer to Ryder fl0|| , for example, for a detailed 
introduction. We decompose the field <f> in terms of a complete set of orthonormal eigenfunctions 
of the vacuum fluctuations 



a n e n (x) + 6 — a n e n (x) 



(43) 



where a is the annihilation and a 1 the creation operator (we neglect the "on them). The hamiltonian 
(0) becomes 

2£f = ^2 Un ( a » a « + a i a n) = X! ujn ( 2a i a « + 1) ( 44 ) 

n n 

using the orthonormality relations of e and the commutation relation between a n and a n . The term 
a n a n is viewed as the number operator N n and gives the number of nth oscillators that are excited. 
We see that the sum J2 n •"■ ^ s divergent and the common procedure is to re-define the hamiltonian. 
We are free to choose the zero of energy and are allowed to neglect the 1. Phrased differently, 
we normal-order the hamiltonian by writing all annihilation operators to the right of all creation 
operators. Thus, we get 

2 : H : = : ^]u; n (a n e4 +a^a„) : = 2^cj„aj i a n (45) 
flip 

where :: stands for normal-ordering. The relations between a normal-ordered and non-ordered 
product of the fields are 

: 4 : = 4 + acj) 2 + f3 

:<j> 2 : = 4> 2 + 5 (46) 

where a, (3 and 7 are constants. We write the normal-ordered hamiltonian as our non-ordered 
hamiltonian plus two counter-terms that arise from the relations 



If : = // - / r/r ( ^6m 2 (j) 2 + ~, ) (47) 
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where Sm is the mass correction to the field and can be evaluated using the one-loop Feynman 
diagram. The constant 7 is not of any importance, because it will cancel itself out due to its 
presence in both, the vacuum and the 1-kink, hamiltonian. The additional term to Efi n i te and the 
divergent term (^) come from subtracting the counter-term of the vacuum from the counter-term 
of the 1-kink and we get 



1 



-5m 



m 

T 



dx 



tanh 



mx 
V2 



V2m 
A 



Sm 2 



(48) 



We evaluate dm 2 by using the equivalent expression in 4 theory. We refer to Ryder (T(], section 
6.4] for a detailed discussion. The standard formula (Ryder: eq 6.95) in perturbation theory for 
the (j> 4 model is 

Sm 2 = ^<?A F (0) (49) 

where g is the coupling and A F (0) the free particle propagator of a loop diagram i.e Ap(x — x). 
We have to be careful when adapting the result to our case. Three modifications to 4> 4 (Ryder eq 
6.65) are necessary: 

• g /4\ = A/4 and g = 6A. 

• The theory should be in (1+1) dimensions. 

• There is only one vacuum. The vacuum eigenvalues are k 2 + m 2 and those of our <fr kink 
theory are k 2 + 2m 2 . Therefore, we need to change the mass m 2 to 2m 2 . 

The (1+1) dimensional free particle loop propagator (Ryder: eq 6.14) with modified mass 2m 2 has 
the form 

. , nS h f dk 2 ,„„, 



(27T) 



2m 2 



where we have a pole at k 2 = 2m 2 and the two- vector k equals (E, —k). We evaluate the double 
integral further 



4tt 2 



A F (0) 



dk 



dE 



E 2 - (fc 2 + 2m 2 ) 
dk 



\/k 2 + 2m 2 



(51) 



where we have integrated over E. We change to the variable p 
Substituting everything into (E9h, we get 



Sm 2 



3A7i 



lim 



dp 



47T A^oo^A ^Jp2 _|_ 2 

which we substitute into the additional term mi 



k/m and put a cut-off A on p. 

(52) 



3.3 Finite Mass Correction 

Finally, we are able to write the quantum mass of the 1-kink as 



M = E K 



E CT — Efi n i te + 



3V2 



Air A- 



lim 



dp 



2(p 2 + 1) 



vV + 2 y/p 2 +2 (2p 2 + 1) 



(53) 
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We have done the integral using a cut-off A with Maple. Both terms produce the same logarithmic 
divergent term which cancel each other out. Taking the cut-off to infinity, we get the final answer 
for the mass of a 1-kink up to first order quantum corrections 



M 



m 

T 



Tim 



1 /3 3 

6V2 ~^7i 



(54) 



where we have written Efi n i te out explicitly. The first term is just the total energy of the classical 
1-kink solution. Note that the presence of 1/ A indicates the non-perturbative nature of the solution. 
To zeroth order in A and first order in h, we have the first quantum correction. It is only valid in the 
weak-coupling limit. The next term of the quantum correction would be of order Xh 2 . Rajaraman 
H section 5.4-5.6] gives a detailed interpretation of the result and also explains why the effect of 
the counter-terms on the kink solution and the zero mode are effects of order A. 

This concludes our derivation of the quantum mass. In the next section, we show that it is 
possible to get a formula for the mass correction which allows us to quantify the contribution of 
the different modes and compute the mass correction numerically by using numerically computed 
lowest eigenmodes. 



4 Trace formula: Derivation 



The trace formula has been first published by Cahill et al. ||, but an explicit derivation has not 
been given in their paper. We derive the formula in this section^. We start by writing out the 
hamiltonian (^) 

27Y(i, x) = ir 2 (t, x) + e(i, x)A 2 i(t, x) (55) 
and the equation of motion of the normal modes is 

e(i, x) = -uj 2 e{t,x) (56) 

with the eigenvalue equation 

A 2 e l (x) = lo 2 6i(x) (57) 

where we ignore the * on the quantum field. A 2 depends on the static solution around which we 
expand. We label Ay the operator for the vacuum and A 2 K the operator for the kink. We expand 
the quantum fluctuation e(t, x) in terms of the normal modes of the vacuum, which we label exit, x), 
and the 1-kink, which we label ty(t, x). In terms of the plane waves of the mesons with eigenvalue 

Wife, 



e v (t,x) = 



ra{k)e lkz 



a\k)e 



—ika 



and, in terms of the normal modes of the 1-kink with eigenvalue uj. 

-iuj n t „iu) n t 



e K (t,x) = Y 



r ^n^n(^') 



(58) 



(59) 



where the plane waves exp(ikx) and the normal modes e n are orthonormal eigenfunctions. The 
next step involves writing the annihilation and creation operators of the eigenmodes in terms of 
the creation and annihilation operators of the planes waves. By definition, 



e K (t,x) 
e K (t,x) 



e v (t,x) 
e v (t,x). 



(60) 
(61) 



3 private communication by Barnes 
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We then integrate over both equations with x, use the fact that the eigenmodes are a complete 
orthonormal set ( ]To| ) and solve the set of two equations for a n and . We obtain 



4 = ^E 

k 



a\k)~e n {-k) 
at(fc)e„(-fc) 



a){k)l n (k) 



UJk 
0J k 



(62) 
(63) 



where e n (k) is the exponential Fourier transform, J dx exp(ikx)e n (x) , of e n (x). We calculate the 
hamiltonian in terms of soliton normal modes annihilation and creation operators in the last section 
(|45| ) and expand the term in terms of the annihilation and creation operators of the vacuum. We 
get 



rt, 



'terms with a^a, aa and aV' 
ic n £a(^ n (fc) (V^~VW a1 



(0?"(-0 



= 'terms with a'a, aa and a'a' 



(64) 



where we have used the commutation relation between operators [a(k),a' (I)) = Sk,i and merged 
the resulting a^a term into the 'collective' term. Finally, we can express the un-ordered term as 
the normal-ordered term i.e. all the terms with a^a, aa and aW and an extra term: 



uj n ala. 



j^e n {k)e n (-k) 



which leads us to the final answer 

: H : 



(65) 



(66) 



= LOna^an + 8m. 

The mass correction 5m can be expressed as a trace over any complete set of orthonormal states, 

(A K - Ayf 



8m 



1 



(67) 



where Ay is the operator of the vacuum and A K is the operator of the kink perturbations. This 
trace formula is finite. Further further discussions see and [jl). 



5 Trace formula: Theoretical Result 

We use the trace formula to calculate the contribution of the lowest discrete modes to the mass 
correction. Cahill et al. only quote the results in their paper [|| We can re- write the trace formula 
( |6^ ) in the following way 

1 f°° 

Sm= --7 / I / dk \r)K,n(x)rjv,k( x )\ 2 K^t 1 - + (68) 

n J -° c 

7 see also M for explicit calculations 
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which reduces for the special case of the zero mode mass correction to 



Stuq 



dk \r]K,o( x )vv,k(x)\ 2 uJk 



(69) 



where r]K,n are the eigenmodes of the kink and rjv,k the eigenmodes of the vacuum. Finding the 
appropriate Fourier transform is the main technical difficulty in solving these kinds of integrals. 
We have used the Maple library inttrans to find Fourier exponential, cos and sin transforms and 
the book on integral tables by Erdelyi et al. §. 

5.1 (f) A kink model 

We have seen that the <^ 4 -kink has two discrete modes (E4l). The zero mode 



3m _ 2 
: cosh 



4V2 



mx\ 

V5J 



(70) 



with lo 2 k o = and the second discrete mode 



VK,l( x ) 



3m . , / mi 



: sinh 



cosh 



mx 

7f 



(71) 



with uj 2 k l — |m 2 are here given in their normalised form i.e. integration of the mode squared over 
x gives one. The normalised eigenmodes of the vacuum fluctuations are 



W,k( x ) 



, ik.r 



/2tt 



(72) 



with eigenvalues ujv,k — k 2 + 2m 2 . Using (B9), we obtain the following integral 



5m 



(...) / dk y/k 2 + 2r 



dx- 



cosh' 



( rax \ 
V2j 



dy- 



,-iky 



cosh' 



( m v \ 

\V2j 



(73) 



which we simplify to 

poo 

Sm = (...)/ dk Vk 2 + 2r 



-1 2 



dx- 



cos(kx) 



cosh' 



( mx \ 
V2j 



(...) 



dk 



\Jk 2 + 2m 2 k 2 



sinh' 



3m 

2V27T 3 



dq 



g 2 v /q 2 + ir 2 
sinh 2 q 



(74) 



/ 



where we have used Euler's formula, the Fourier cos transform of cosh 2 and changed to the variable 
q = -j^- Using (68), we obtain the correction to the mass from the second discrete mode 



Smi 



(...) 



dk 



3m' 



(k 2 + 2m 2 )-i - 3m 2 + yjk 2 + 2m 2 



dx sin(fca;) 



sinh(^| 



cosh' 



( mx \ 
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h 

where we have used Euler's formula, the Fourier sin transform of sinhcosh -2 and changed to the 
variable q = -jS—. 

We have evaluated Iq and I\ numerically with Maple and obtain 

h = 11.247 

-h = 0.827 
4 

and the contributions to the mass correction are 

6m + Smi = (-0.384- 0.0283)m = -0.413m (76) 

compared to the full quantum correction (pij) 

8m = -0.471m (77) 

where we have set Ti to one. Finally, we find that the zero mode contributes 81.5% and the second 
discrete mode 6%, thus in total 87.5%, to the total mass correction. 

We are also interested to what value of k the normal modes of the vacuum fluctuations have to 
go to give a reliable answer to the mass correction. Intuitively, the zero and second discrete mode 
are localised and its norm with the long wavelength vacuum modes should become very small. We 
have put a cut-off A on our integral Iq and jli and evaluated the integrals as a function of the 
cut-off. We have done this numerically with Maple. Figure |l| shows that we only need to go up to 
a cut-off of around q = 5 for the zero mode and of around q — 15 for the second discrete mode. 
(k= y^HLq) This is good news, for we can ignore long wavelength vacuum modes. 



5.2 Sine-Gordon model 



We do the same for the Sine-Gordon model. There is only one discrete mode, the zero mode of 
translation: see section ^ and Q . The normalised zero mode has the form 



Vk,o(x) = y^cosh X {mx) 



(78) 



with eigenvalue lok.q = 0. The eigenmodes of the vacuum fluctuations are the same as before, but 
the eigenvalues change to u>v,k — k 2 + m 2 . The mass correction takes the form 



5mo = (...) y dk\J k 2 + m 2 

= (...)/*v^co S h 



dx cos(kx) cosh 1 (mx) 



nk 
2m 



cosh q 



(79) 
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Figure 1: Value of Iq (upper curve) and I\ (lower curve) as a function of their cut-off A for q 



where we have used Euler's formula, the Fourier cos transform of cosh 1 and changed to the variable 
a = — 

We have evaluated the integral numerically with Maple and obtain 



Tfl Tfl 

Sm = - 3.572— = - 0.893 — = - 0.284 m 

4-7T 7T 



compared to the full quantum correction (||, |J) 

Tfl 

6m — = — 0.318 m 

TT 



(80) 



(81) 



where we have set h to one. Finally, we find that the zero mode contributes 89.3% to the total 
mass correction. 

Again, we are interested to what value of k the normal modes of the vacuum fluctuations have 
to go to give a reliable answer to the mass correction. We have put a cut-off A on the integral 
and evaluated the integral as a function of the cut-off. We have done this numerically with Maple. 
Figure || shows that we only need to go up to a cut-off of around q = 5 for the zero mode. We do 
not need to include long wavelengths vacuum modes. 



6 Trace formula: Numerical Result 

In the last section, we have applied the trace formula to the <f> kink and the Sine-Gordon model. The 
results are clear-cut. The contribution from the discrete modes to the quantum mass correction 
is dominant (more than 80%). Further, we do not need to probe our discrete modes for long 
wavelength of the vacuum mode. This is good news for numerical methods and we can limit 
ourselves to the lowest normal modes of fluctuations in both the vacuum and the kink sector. 
Moreover, in (1+1) dimension, we are not restricted by memory or computational needs and can 
include all the vacuum and kink modes. 
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Figure 2: The value of the integral versus its cut-off A for q 



6.1 Preparation 

We calculate the mass correction for the 4> 4 kink and Sine-Gordon model. We have set the mass 
m = 1 and coupling A = 1 for simplicity. The energy functional has the form 



for the Sine-Gordon model; where (f)' = 4j and <^> = Using appropriate boundary conditions, we 
find the minimal-energy configuration, which we call <fi stl for both models in the topological charge 
sector one. We use the Gauss-Seidel overrelaxation method. Our box size is L = 40 from -20 to 20 
and we use 1600 points. Thus, the lattice spacing is dx = 0.025. 

The corresponding eigenvalue vacuum and kink operators in terms of the static solution (see ^) 




(82) 



for the c/> 4 kink and 




(83) 



are 




(84) 



for the (jft kink and 



4 2 




K 




(85) 



14 




for the Sine-Gordon model. In Numerical Techniques, we describe three different methods that 
solve the discrete eigenvalue problem. The trivial matrix diagonalisation is the more accurate and 
simplest one. However, we have to admit that the computational time grows as the cube of the 
number of points and the technique cannot be used in two or more dimensions. We substitute the 
value of the numerically minimised field^ into the discretised eigenvalue equation and diagonalise 
the resulting matrix with periodic boundary conditions. 

Figure shows the first four normal modes of the 4> 4 kink. The numerical eigenvalues are 
—5.78 10~ , 1.49998, 2 and 2.03072 compared to the exact eigenvalues (zero mode), 1.5 (first 
mode) and 2 (start of 'continuous' modes). Figure |] shows the first four normal modes of the Sine- 
Gordon soliton. The numerical eigenvalues are 5.12 10~ 8 , 1.00682, 1.00682 and 1.06128 compared 
to the exact eigenvalues (zero mode) and 1 (start of 'continuous' modes). 

6.2 Results 

We have all the information needed to compute the mass correction. The trace formula has the 
form 

n,k \ i / 

which is similar to d6§|). rjK,n( x i) refers to the nth eigenmode of the kink, rfv,k{%i) refers to the kth 
eigenmode of the vacuum and Xi to the position of the ith lattice point. 

We start with the 4 kink. We sum up the contributions to the mass correction, mode by 
mode. Figure ^| shows that the mass correction approaches an asymptotic limit. The first few 
mode contributions are the most important ones. We only include mode contributions up to mode 

we used the relaxation method described in Numerical Methods 
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Figure 4: The Sine-Gordon eigenmodes from zero to three 
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Figure 5: The total contribution to the mass correction mode by mode 
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Figure 6: Norm between four kink modes and the vacuum modes 



200. Figure |6j shows that the norm is mostly zero and peaks for the norm between the nth kink 
mode and the nth vacuum mode. For the lowest kink mode contribution, only the lowest vacuum 
modes are important. For the highest kink mode contributions, only the highest vacuum modes 
are important. The discrete mode corrections are 



Smo — 
Smi — 

which are very close to their exact values fl76|). 

8m — 



-0.384626 

-0.0282964 (87) 



- 0.471097 (88) 



The numerical value of the mass correction is 



compared to the exact value of -0.471113 and is 99.997% accurate. This is a very satisfactory result. 

We turn our attention to the Sine-Gordon kink. We sum up the contributions to the mass 
correction mode by mode. Figure |?] shows that the mass correction approaches an asymptotic limit. 
The first few mode contributions are the most important ones. Figure |^ shows the contribution 
of the first mode for each vacuum mode. Note that some contributions are zero, because the first 
mode is an odd function and some of the vacuum modes are even functions. The discrete mode i.e. 
zero mode correction of the Sine-Gordon kink is 



6m = -0.28402 (89) 
which is very close to the exact value ( |S0"| ) . The numerical value of the mass correction is 

5m = -0.318144 (90) 
compared to the exact value of -0.318309 and is 99.95% accurate. This is a very satisfactory result. 
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Figure 7: The total contribution to the mass correction mode by mode 
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Figure 8: Mass correction of the first kink mode versus the vacuum modes 
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7 Conclusion 



We have shown that the trace formula works very well in (1+1) dimensions for the Sine-Gordon 
and the <j) A kink model. The numerical quantum mass correction is very close to the exact one. Our 
technique can be applied with ease to any (1+1) dimensional theory. This allows us to calculate 
the mass correction to non-integrable solitonic systems, for example. Specifically, we are interested 
in the Sine-Skyrme model ^ and plan to study the mass correction numerically. Or, we can look 
at the mass correction of multi-solitons, for example. 

There are two drawbacks. We have used a brute force matrix diagonalisation to find the eigen- 
values. It works very well and is reasonably fast in (1+1) dimensions. However, if you are to 
implement the trace formula in higher dimensions, you will have to use a different technique. 
Barnes and Turok in 0] have used a diffusion equation method; as discussed in the chapter on 
numerical methods. They compute the zero mode, then project it out of the initial configuration, 
get the next mode and so on. Only the first modes are accurate, because the errors are summing 
up. Computational restrictions also limit the calculations to the first few modes. 
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